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SDMUAEI 


A metkod. is presented for calculating tke aerodynamic loading, tke 
diTergence speed, and certain stability derivatiTes of vings and tail 
surfaces of arbitrary plan form, and stiffness. Provision is made for 
using either stiffness curves and root-rotation constants or influence 
coefficients in the analysis. Collating fozms, tables of numerical 
constants rbquired in the analysis, and, an iUuHtratlve ezas^le are 
Included to facilitate calculations by means of the method. 


USTROEOCTION 


Ibe distribution of the aerodynamic loading on vings and tall 
surfaces is Important both for the structural analysis of these com— 
ponents, since it determines the applied bending moment and torque 
acting at any station, and for their aerodynamic analysis, since it 
affects the stability derivatives to a large extent. At hi^ speeds the 
aerodynamic loading, particularly in the case of svept vings, is 
greatly affected by the structural deformations caused by the loading. 
33ie present paper is concerned vlth the determination of the effects of 
structural flexibility on the aerodynamic loading of wings of arbitrary 
plan form and stiffness. 

The problem of load distribution was analyzed for imswept flexible 
wings as early as 1926 (refereiice l) but has received relatively little 
attention since. The only new effect considered in subsequent work is 
aerodynamic Induction (reference 2), No work appears to have been done 
on the loading of flexible swept wings. The related problem of aero- 
elastic divergence of swept wings with certain prescribed stiffness 
variations has been treated In reference 3* 

Eie present paper ti*eats the problem of aerodynamic loading by 
matrix methods. Aerodynamic induction is taken into account approxi- 
mately, since sixitable aerodynamic influence coefficients are not avail- 
able for wings of arbitrary plan foim. Wh«i they become available they 
can readily be Incorporated in this method. Structural flexibility is 
taken into account in the form of either calculated stiffness variations 
or measured influence coefficients. The required integrating matrices 
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Etre presented for totli a six— point and a ten-point solution. For the 
six— point solution convenient computing forms are included, as well. The 
method is illustrated 137 means of an example. In ad.dltlcsn to the anal 7 sls 
of the aerodynamic loading, the determination of the related divergence 
speed and of certain stahility derivatives is discussed. 

For the convenience of the reader unfamiliar with matrix terminology 
a summary of matrix methods has "been included in the appendix. The 
sections on "APPUCATION OF TTTC METHOD" and, in particular, ’'Instructions 
for Solution" may he read without reference to the section "DjJitdLVATIQN OF 
THE METHOD." 


SYMBOIS 


A 

[a] 


a* 

ac 

b 

c 

c 

% 

“Hr 


aspect ratio 
aeroelastic matrix 


cos^As^2e^ Cj.2'' 

dimensionless parameter i ^ 


(GJ). 


parameter cos As^^p) 


section aerodynamic center, measured from leading edge, 
fraction of chord 

wing span. Inches 

chord measiared parallel to the air stream, inches 
average wing chord, inches 

section lift coefficient ^ 

\<ic/ 

/W 

wing lift coefficient — 

\<lS, 

rate of change of damping-moment coefficient in roll with tip 
helix angle 

/M* \ 

wing root hending-moment coefficient ( — 

\qSb/ 

wing rolling-moment coefficient 
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d 


dimensionless parameter 


mgg. cos^s^3oj 


tan A 


(EE)-, 


■) 


EH tending stiffness In planes perpend/ oular to the elastic 

axis, pound— Inches^ 

ea location of elastic axis measured from leading edge, fraction 

of chord 


e^^ distance from reference axis to section aerodynamic center 

(positive forward) (see fig. l), fraction of chord 


W 

w 

[%] 

w 

K] 

H 

N 

N 


torsional stiffness in planes perpendicular to the elastic 
axis, pound— Inches^ 

unit matrix 

matrix defined hy equation (12.) 

Integrating matrix for sin^e integration from tip to root 
first row of Kj_ matrix 

Integrating matrix for double integration from tip to root 
first row of malarlx 

integrating matrix for single integration from root to tip 

matrix relating concentrated and accumulated torque 

matrix relating concentrated loads and accumulated bending 
moments 

matrix converting torques due to distributed loads to torques 
due to concentrated torques 

matrix converting bending moments due to distributed loads to 
bending moments due to concentrated loads 


L lift, pounds 

lift on both wings but excluding lift on part of wing covered 
by fuselage, pounds 

I running air load along the reference axis, pounds per inch 

M accumulated bending moment (in planes pez-pendloular to the 

reference axis imless specified otherwise), inch— pounds 
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effective section lift— curve slope for angles of attaclc due to 
deformation, per radian 

eIq^ effective section lift-curve elope for additional —type angle— 

of-attack distributions, per radian 




P 


Qr 

1 


E 


s 


®A 

T 


w 


7 


section liffc-cuive slope, per radian 
concentrated load, pounds 
root— twist constants (see equation (9)) 
root-bending constant (see equation (9)) 
dynamic pressure, pounds per square inch 
concentrated torque, Inchr-pounds 

total wing area including part of wing covered by fuselage, 
square inches 

distance from wing root to wing tip peiT>endicular to the air 
stream (see fig. l). Inches 

length of wing along reference axis (see fig. l), inches 

accimilated torque (in planes perpendicular to the reference 
axis unless specified otheivise), inch-pounds 

distance between the effective root aiid the innermost conqjlete 
section of the torsion box perpendicular to the elastic 
axis. Inches 

lateral ordinate measured from wing root. Inches 


a 


a 


angle of attack, radians 


equivalent angle of attack, radians 



r local dihedral angle due to deformation or slope of wing 

deflection curve at reference axis, radians 

^ structural deflection, inches 

T) distance along reference axis, inches 

angle of sweepback (measiired to the reference axis unless 
specified otherwise), degrees 


A 
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q> 


influence ooeffioient matrix for wing twist in planes parallel 
to the air stream due to concentrated unit loads applied at 
the reference aii.s, radians per pound 

influence coefficient matrix for wing twist in planes parallel 
to the air stream due to concentrated unit torques applied 
in planes parallel to the air stream^ radians per inch- 
pound 

angle of twist in planes perpendicular to the reference axis, 
radians 


Subscripts : 
c midchord 

D divergence 

fw flexible wing 

g geometric 

LB leading edge 

M due to bending moment 

MA.C pertaining to the mean aerodynamic chord 

r at root or effective root 

rw rigid wing 

s structural (due to structural deformations) 

T due to torque 

H!B trailing edge 


w wing exclusive of fuselage 

Prime mark: 

• in or pertaining to seoticms parallel to air stream rather 

than perpendicular to the reference axis 


Matrix notation: 


{} 

LJ 


column matrix 


row matrix 
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[] 

square matrix 

[°] 

diagonal matrix 


nEEITAHOIf OF THE METHOD 
Method En 5 )loylng Stiffness Ouarves 


Assungrfcl OUR In the development of the method the following 
assun^tions are made: 

(a) The effects of eierodynamlc induction may be taJcen into account 
by applying an over^-ell correction to strip theory and rounding off the 
resulting load distribution at the tip. 

(b) All deflections and angles of attack are small. 

(c) The wing is mounted fleslbly at an effective root peip)endlc-ular 
to the elastic axis throu^ the Intersection of the elastic axis and 
the fuselage (see fig. l), the root rotations being proportional to the 
root bending moment and root torque, 

(d) An elastic axis exists in the outer portion of the vlng^ this 
axis being defined as the elastic axis the wing would have if it were 
mounted rigidly some distance outboai*d of the root approximately perpen- 
dicular to the midchord line. (Hear the root the elastic axis is defined 
as the extension of the outboard elastic axis.) 

(e) All deformations are given by the elementary theories of 
bending and of torsion about the reference axis, which in this case is 
the elastic axis. 

Air loads .— In keeping with assumptions (a) and (b) the force on 
a wing section of unit width parallel to the direction of fli^t is 

7* = q.c cos A(^oce 4 m^a^ 

= meq.c cos Aa (l) 


where the equivalent angle of attack a is defined by 


_ “®1 


(la) 
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The effective seotion lift— omrve slope tOq tar angle— of-attaolc distribu- 
tions due to structural defozmatlons has been given in reference 3 OQ- 
the basis of the reasoning of reference for subsonic speeds as 




A 

A + 4 cos A 


( 2 ) 


Similarly, the effective lift— curve slope for addition a l — type 

angle— of-attaoh distributions is determined by tl^ same reasoning as 


inei 


= mo 


A 

A. + 2 cos A 


(3) 


Both slopes must be multiplied by cos A, as in eguation (l), in order 
to apply to loads and angles of attack acting on sections and measured 
in planes parallel to the direction of f1 Ight. 

The torque of the running load 2* about the reference axis is 2*e2_c 
for uncambered sections (for cambered sectioaos the torque at zero lift 
must be added and the analysis of the foUovlng paragraphs modified 
accordingly) . Biis torque may be resolved into a running torque about 
i^e elastic axis and a running bending momait about a line perpendicular 
to the elastic axis. HSie running load, torque, and moment must then be 
•multiplied by cos A to yield •fcheir values per unit length along the 
elastic axis, so that 


2 = m©q.o oos%5 (^l-) 

or, in matrix notation, 

= m^q. cos?A[cj (^) 


3!tie running torque and the running bending momiant aire, respectl-vely, 
^20^0 cos A^ and 2e2_c sin a|. The running bCTtdIng mcanent leads to 

accumulated bendizg momehts -Hhich have to be added to *the accumulated 
bending 'moment due to the zunning load. 

The accumulated torque T is obtained from the running torque by 
an integration Inboard from the tip. This integration may be performed 
by a matrix Dhich is based on Siitqpson*B znile •wl'th a modification 

su gg ested by Y. M. Ealkner at the tip. (See append It. ) The effect of 
Jalkner's modification is to roui^ off •the calculated load distribution 
and cause it to go to zero •vrl'th an infinite slope at the tip, as *the 
aerodynamic lift distributions actually do. The matrix is given ^ 
table 1. 
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S1iti11arl7, tlie accxnoulated 'bending moment Is o'btalned a double 
Integration inboard fj^cm tbe tip of tbe running loed. and a sln^e 
integration of the running moment. The double integration may be 
performed by another matrix [Kg] (given in table II ), TJhich is based 

on the equlTalent of SlTi^Bon*s rule for moments, Falkner*s modification 
again being made at the tip. Ihe derivation of the integrating matrices 
is discussed in somewhat greater detail in the appendix. 

The accumulated torque and bending moment may then be written as 

W = ®a[%] {^®1° 

O 

{s} (5) 

« 



W = W - °a[%] {»«!■= a] 



Equations of equilibrium .— Bie equations of equilibrium, of a 
deformed wing referred to the elastic axis are 


GJ 



(7) 


El — = M (8) 

dy 

These equations must be Integrated outboard from, the root to obtain (p 
and F. niie Integratians may be performed by a matrix [Kgj (see 

table HI and appendix), also based on Simpson’s rule without the tip 
modification, however, since the torques and moments go to zero with 
finite and zero slopes, respectively. To the deformations obtained in 
this manner the rotations due to the root deflection, 9^. F^^ 
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mufffc te added. TLe rool; rotatlcais etre defined b7 four dlmsaislonleBS 
constants : 






I- 


v/(GJ) 


(9a) 




w/(Gj)j. 


(9b) 




Qp — 

v/(El)^ 


(9c) 


*^r, 




M 


v/(El)^ 

■wiilch. may be ccaobined Into two other oanstants 


^CtiJ 


W' 


T_ 


v/(GJ), 


= ^ - 


(GJ)r ) 

- — -i- tan AQr. 
(EE)^ ^tJ 


cos A 


(91) 


(9e) 


^ v/(GJ)^ 




(GJ)r A 

— tan AQiJ oos A 


(9f) 


w being defined as In figure 1. !Ihe defoHnatlonB may thea be written as 



O 


1 ba 


W * f «%0k]{“} 


( 10 ) 
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“ 

- 

0 

- 


II 


N 

'(EE)r' 

EE 

- ^ 


r 

_ 


L. J 



vhere the matrix 



is defined by 



( 11 ) 



0 0 0 0 
10 0 0 
10 0 0 
10 0 0 
10 0 0 


( 12 ) 


TOie angle of attack due to the structural deformations Og Is 
related to qp and P 


ctg = (<p — P tan A) cos A 


(13) 


If equatlans (5), (6), (lO), and (U) are substituted in the matrix 
equivalent of equaticsn ( 13 ), the following relation is obtained 


{®s} “ 

where the aeroelantlc matrix [aJ is defined by 


(14) 


W = 


O 


ig; [Igi.] 




aN [tT-J 


W 


^A 


Sj, e-) ^ Cr cos A 
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and 'Uie parameters 


m^q. oos^s^2ei^Cj.2 

(GJ)r 


(16) 


and 


d 

a 


(Et)j, oos A 


tan A 


(l6a) 


are substantial I7 the same parameters as those -used In reference (3). 

Solution of the eguatlgns .— If It Is desired to calculate the 
aerodynamic loading corresponding to a glren geometrical an^e— of— attack 
distribution and dynamic pressurej equation (ih) may be refwrltten as 
f oUovrs : 



(IT) 


In this form It constitutes a set of linear simultaneous equations for 
the a Talues In terms of ctg Talues. Stom the ceilculated a values 

the lift distribution may be determined from either equation (l) or (i»-). 


Ihe divergence dynamic pressure may be obtained .from, equation (I7 ) 
by setting the determinant of the square matrlz on the left side of the 
equation equal to zero. iEhls procedure Is equivalent to setting Og 

equal to zero In the term a of equation (l 4 ), so that 


|oaj> = a[A]|oB| (18) 

!Ehe critical value of a Is then determined by matrix Iteration and 
hence the divergence dynamic pressure frcm equation (I6). 


Method En^jloying Influence Coefficients 

The assunq>tlons made in the preceding sections concerning the 
behavior of the vlng structure are unnecessary if influence coefficients 
for the glvKi structure are available from, test data or refined methods 
of calculation. The coefficients most convenient for this analysis are 
those giving the rotation of the structure in planes parallel to the 
directicai of f 1 1.^t due to vertical loads applied along a convenient 
reference axis and due to torques about lines perpendicular to the 
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direction, of fll^t. Since It is usioally more coavenlent to appl 7 
concentrated rather than dletrihuted loads In Btructiiral tests, the 
Influence coefficients vlll he considered in this analysis to have 
heen obtained in this manner. 

The angle of structural deformation ttg may he eaqjressed in terms 
of the influence coefficients ^ and 4^ as follows: 

{“»] - MW * [*p]W (w) 


where the E*s and P*s are arbitrary concentrated torques and loads, 
the latter being applied at the reference axis. Ihe accumulated torques 
and bending momiints about lines pe^pendlc^llar and parallel, respectively, 
to the direction of fli^t may be related to the concentrated torques 
and loads by means of the surnmation matrices 

appendix) as follows: 


{t.} . [E4.](e} - a{m.} 

{“■} - 


( 20 ) 


( 21 ) 


These relations may be solved for the vsG.ues of E and P required to 
produce given distributions of accumulated torque and bending moment 



W = 


( 22 ) 

( 23 ) 


The accuamlated 
load are then 


torques and bending moments produced by the air 





tan A 


( 24 ) 



( 25 ) 
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TTpcoi substituting eqmtlcms (22), ( 23 ), ( 25 ) and (l) Into 

equation (19 "tlie following equation Is obtained 

|ag^ = a-*[*-*] (26) 

•wiiere 


a* = m© cos AqSyOy ( 27 ) 


M 




(as) 


wiiere. In turn 


[%] - w "■[%] 

( 29 a) 

_ ~ 1 . .. 

M = f5j H 

( 29 b) 


are given In tables lY and V. 

Ihe solution of equation ( 26 ) Is obtained in the manner previously 
described for equation (ill-). 


APPLICAnOir OB* TUT?; ME33HOD 
Pet anal nation of the Structural Parameters 


At the time an aeroelastic analysis is performed no experimental 
stiffness data are usually available, so that the calculated stiffness 
curves must be used. In carder to use Idieee curves it Is necessary to 
assume the existence of a reasonably stral^t elastic ails. The location 
of this axis may be estimated by considering It to be the line connecting 
the ^ear centers of the individual sections. If the elastic axis 
obtained In this manner 1s not reasonably stral^t within a few percent 
of the chord, the results of the analysis may not be sufficiently reliable. 
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Hae stiffnesses GJ and El do not hare muoli physical, significance 
inboard of the last point where there is a complete cross section of the 
torsion box. (See fig. 1. ) In order to arrive at estimates of the root 
stiffnesses (OJ)^ and (iu)p, which serve primarily as reference values 

in this analysis, the stiffness curves have to be extended. It is con- 
venient to consider the stiffnesses to be constant inboard of the last 
con 5 )lete section of the torsion box; this procedure should yield conservar- 
tive values of the root rotations. 

Ihe most difficult problem, incurred in analyzing the deflections on 
the basis of stiffness curves appears to be the estimation of the root 
rotations. As used in this analysis, they are the torsion and bending 
deflections ln^josed by the triangular inner portion of the wing and the 
carry— throu^ bay on the rest of the wing. As seen in figure 2, which 
is plotted from the data of reference 5, these values are essentially 
constant along the span, so that they actually constitute rigid- body 
rotations. (The bending rotations have been obtained by taking the 
difference in slope between curves calculated by considering the wing 
to be cantilevered at the effective root — the root used to calculate 
torsional deformations in reference 5 — and the averages of the leading- 
edge and t railing-edge deflections actually meaieured. The twists were 
obtained by subtracting the twists calculated on the basis of the assumed 
effective root from the measured twists. ) 

The rotations should in any practical case be calculated by analyzing 
the triangular root and the carry— through bay and made dimensionless by 
means of equation (9). If such an analysis is not available, the dimension- 
less rotation parameters shown in figure 2 may be used as a guldej it 
must be kept in mind, however, that in the case of a sweptforward wing 
the parameters Qq^ and Qjjj, would have the opposite sign. 

Once the stmcture under investigation is bullt^ fairly slE5>le 
deflection tests, similar to those performed in reference 5j “*7 "be used 
to check the root— rotation parameters by calculating the differences 
between the observed rotations and those calciilated by simple beam 
theory considering the wing cantilevered at the effective root; at the 
same time the existence and estimated location of the elastic axis may 
be verified. If the eaperimental program is fairly extensive it is 
desirable to measure influence coefficients directly. These Influence 
coefficients can then be used in conjunction with the alternate method 
described In the preceding section to obtain a quick check on the aero- 
elastic analysis based on calculated stiffnesses. 

The influence coefficients used in the analysis consist of the 
rotations of sections parallel to the direction of fli^t due to 
concentrated unit torques in planes para.llel to the plane of symmetry 
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or conoentrated ■unit loads at tlie reference line, 
radians are entered in tatles of the form; 


[*b] 

TWIST AT STATION 7 /s^ IJOE TO 
UNIT CONCENTRATED TORQDE AT 


7i/s^ 


7 i/sw 

y/Bw 

0.1 

0.3 

0.5 

0.7 

0.9 

1.0 

0 






■ 

0.1 






■ 

0.3 

■ 

■ 

■ 

■ 


■ 

0.5 

■ 

■ 

■ 

■ 


■ 

0.7 

■ 

■ 

■ 

■ 


■ 

0.9 

■ 

■ 

■ 

■ 


■ 


These rotations in 

[•p] 

TWISP AT SEATION 7 /&W HIE TO 
UNIT CONCENTRATED LOAD AT 


7l/®w 


Bi 

0.1 

0.3 

0.5 

0.7 

0.9 

1.0 

0 







0.1 







0.3 







0.5 







0.7 







0.9 








These particular tables would he used for a six— point analysis; similar 
tables would be used for a ten-point analysis. In either case it is 
to be noted that the twists are measured at values of y/s^ from. 0 to 0.9^ 

whereas the loads are ai^ilied at yi/s^ values from 0.1 to 1.0. The 

tables obtained in this manner constitute the desired influence- 
coefficient matrices. 

If the wing sections are found to twist nonunif oimLy, so that they 
become cambered in effect, the angles of twist Og to be entered in 
the influence-coefficient matrices have to be defined in a different 
manner according to whether the aeroelastic analysis is perfozmed for 
subsonic or supersoonic speeds. At subsonic speeds the lift depends cn 
the elope of the mean camber line at the three-quarter-chord point, so 
that the effective tngle of attack is 


Og = 


(^c 


c 


( 30 ) 
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At supersonic speeds the lift depends primarily on the avera^ slope of 
the mean camber llne^ so that 


^5CE 


a-, = 


(31) 


Determination of the Aerodynamic Parameters 

The selection of the aerodynamic parameters m^ and e^ for the 
calctilatlon of the divergence speed has been disctissed in reference 3« 

For ceQ-Culating the aerodynamic loadli^ at a given fli^t condition the 
aerodynamic parameters are chosen for that fll^t condition. The 
effective lift— curve slopes biq and nig^ are applicable only to subsonic 

subcritical speeds. At higher speeds no simple span correction is 
available; neglect of the span correction tends to be conservative for 
calculation of the divergence speed and the aerodynamic loading, however. 


Instructions for Solution 


Two sets of integrating matrices have been prepared for a eli-point 
solution and a ten— point solution, respectively. The former should be 
adequate for mi practical purposes; only where the stiffness curves are 
very Irregular near the root does the ten— point solution have to be 
resorted to. The points considered by the two sets of tables are at 



0.2, 0.4, 0.6, 0.8, and 0.9 for the shorter solution and 
0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9 for the longer 


solution. The procedure to be followed for either solution is id^tlcal; 
althou^ computing forms are presented in this paper only for the six- 
point solution, their extension to apply to the ten- point solution Is 
obvious. 


Calculation of the matrices .- The first step in the aeroelastlo 
analysis by means of the stiffness curves is the calculation of the 
aeroelastic matrix [aJ from the physical and geometrical parameters 
of the wing. These parameters are conveniently tabulated in a foim of 
the type shown in table VlCa). The computation is then carried out 
according to the instructions of table Vl(b), each step in the procedure 
being identified by the. number in the upper left comer of each box. 

It must be kept in mind that many of the operations call for matrix 
multiplications where the order of the multiplicands is of Inqportance. 

(a brief summary of matrix methods is presented in the appendix. ) The 
aeroelastic matrix is obtained as the last step (step 13 ) of the 
computations in this form. 
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A special case arises ■wiien e-j_ is zero. 


If 


is not zero 


along th.6 remainder of the span^ its value at^some point other tTinn the 


root may be used as a reference value. The 


fl 

e 


Ir 




:) 


matrix, the 


multiplying factors of steps 8 and 9 as well as the definition of the 
parameter a are then based on this other reference value rather 
than 01 • If ©n 


.-r - ®1 

to 8 may be emitted. 


is zero along the entire span, step 1 and steps 3 
Steps 9 to 13 should be modified as follows; 


Step 9 

Step 10 

Step 11 
Step 12 

Step 13 


w. Si [r 1 

(GJ) s^tanA*-®-* 

r 

[®]- [®] 

As is 
Qmlt 

We^.0 = [®] [®] 


If influence coefficients of the proper type are available, the 
calculation of the aeroelastic matrix [A*J Is carried out directly 
by means of equation ( 28 ). 

Solution for divergence dynamic pressure .— In order to determine 
the -roluft of the pan ejne ter a or i7 corresponding to divergence, the 
aeroelastic matrix [a] or [a*] Is Iterated (see appendix). Iteble YH(a) 
may be used for this purpose. Hie result is the critical value of a. 
or a* . Hie divergence dynamic pressure is then calculated from, equa- 
tion (l6) or ( 27 )* It is to be noted that this pressure will be in 
pounds per square inch. Since the aeroelastic matrix is Independent of 
the Mach number, except insofar as e^^ varies with Mach number, the 

same critical value of a may be used to calculate the divergence dynamic 
pressure for an entire range of Mach numbers. If the value of ©j_ 

changes, however, as it does between the subsonic and supersonic region, 
the critical value of a has to be calculated for both values of e^* 

If the value of ej_ is zero along the entire span and the [A] 

matrix has been calculated according to the modified I ns tructions, 
iteration of the matrix will give the value of the parameter d at 
divergence. IProw. the definition of d the divergence dynamic pressure 
may then be calculated. 
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Solutlcoa for aerodynamic loading .— In order to calculate the aero- 
dynamic loading corresponding to a given rii^t condition and geometric 
angle— of-ettack distrihution the aeroelastlc matrix [aJ or (A*J is 
multiplied by the value of a or a* csilcula t ed for the given fli^t 
condition and subtracted from the unit matrix £l3* The result may be 
entered in table VH(b). Again it must be noted that the value of the 
aeroelastlc matrix varies with the fll^t condition if S 2 _ varies, so 
that the aeroelastlc matrix corresponding to the proper e^^ value 
must be selected. The resulting matrix constitutes the coefficients of 
a set of simultaneous linear algebraic equations for the unknown values 
of the effective angle— of-attack distribution of the deformed wing 

1^^— in terms of the known angle— of-attack values of the rigid 

wing Table Yll(b) is set up for the ceilculation of the additional 

loading, the damplng-in-roll loading, and a third etrbitrary loading; as 
many loadings as desired may, of course, be calculated by this method. 

The solution of the equations may be carried out in any convenient manner. 
The form of table Vixfb) has been prepai^ed for use in conjunction with 
Crout*a method of solving linear simultaneous equations (reference 6). 

In the case where en is zero along the span, the headings at the 
top of table VTl(b) should be modified to read 



p] - 0 

L J 


viiere has heen calculated according to the modified instructions 

and d has”’been obtained by iterating 


The values of 


fe-} 


calculated for the additional load 
constitute values of the ratio 




distribution (ag = l) 

forces. The section loading of the flexible wing is obtained from the 
relation 


fw/*^Vw 

in view of the assungstlons made concerning the air 


“I “ y 


(32) 
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or in dimensionless form. 

£21 = m _2-/^ 

°r ®1 Cr\iBej^ ) 


die wing lift coefficient defined by tbe relation. 
anri the wlng bendlng-mcment coefficient defined by 


( 32 a) 


(33) 


9bm,- 


“!e 

qSb 




may be obtained by integrating the load distribution. These Integrations 
may be performed conTeniently by multiplying the co^/cp Talues by the 

first rows of the fEjl and [KgJ matrices, respectlTely. Thus 





(35) 


and 



(35a) 


SNilsr} 


(36) 


•=H>ei 



(36a) 
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The lateral center of presstire of tlie wing load 
determined from, the relation 


— may then be 




( 37 ^ 


The fore-end— eift location of the aerodynamic center of the wing load 
measured rearward of the leading edge of the mean aerodynamic ^ <^hord as 
a fracticxn of the mean aerodynamic chord may be estimated frcan the 
relation 


(ac 


M.A.C. 


ac + 


y-y-MAc 

M.A.C. 


tan 


(38) 


where Is the sweep of the section aerodynamic center line. 

For any other geometrical angle— of -attack distributions such as 
those due to built-in twist or those due to rolling, the aanie section 
lift— curve slope shoitld be used as for the structural deformations, so 

TITIq .- 

that is replaced by mg and ^ — is unity in eqTiatlons (32) to (38). 

For the damping-in— roll distribution with a tip helix angle of 1 radian 


‘"g = 



(39) 


The rolling-moment coefficient due to the wing load is defined by 


C 


2M* 

= fL-n 

Sr qSb 


(1^0) 


It Is seen to be twice the wing bendlng-^noment coefficient. 

The contribution of the wing to other stability derivatives may 
be obtained similarly by integrating the load distributions due to the 
angle— of -attack distributions caused by the motion under consideration, 
as deociribed in reference 7; in the case of swept wings, particular csire 
must be taken in selecting the proper angle— of -attack distribution and 
in accounting for the lateral inclination of the lift vector (see 
reference 4). 

If the aerodynamic loading or the stability derivatives are to be 
obtained for a wide variety of fll^t conditions, it is convenient to 
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S7stematlze the caloulatlons in the follow±ag manner. The aeroelastlc 
matrix is coii5)uted tar "both the sahsonic and supersonic aerodynamic 
center values and iterated for "both cases to obtain the suisonic and 
supersonic values of the divergence parameter ajj. Erom these values 

the divergence dynamic pressure may he ccmputed hy means of equation (15) 
and plotted against Mach number, as su gg ested in reference 3j 0 °. the 
same plot values of the actual dynamic pressure may he plotted against 
Mach number for various altitudes of interest. Such a plot for a ving, 
the physical characteristics of which are given in figure 3, is shown in 
figure it-. 

Since at a given Mach number the ratio a/aj) is equal to the 
ratio q./q.j), the range of h/a^ values of interest may be established • 

from, this plot for both the subsonic and the supersonic region. Several 
representative a/a^ values may then be chosen within the given ranges 
and the corresponding values of a computed from, the previously calculated 
ajj values. The aerodynamic loading is calculated for these values of a 
using the appropriate IaJ matrix and plotted in the form of 
(oc-jT- 3 wibh the ratio a/a^ as a parameter. I*roim these 

curves or from the a values the lift coefficients may be obtained 

/ 

in the form and plotted against a/ajj or 'fclie 

other coefficients may be obtained and plotted in a similar form. 

For any specific fli^t condition the valxie of a/ajj may then 
be obtained from the plot of g and qjj against Mach number. The 

loading, lift coefficient, or other item of Interest m^y be obtained 
from the plots which give these items in terms of the rigid-wing values. 
Once the rigid-wing values at the given Mach number are known, the 
flexible-wing values may then be obtained Immediately. 


Illustrative Exanqple 

In order to illustrate the method described in the preceding 
sections, a typical swept wing has been analyzed. The physical 
geometrical parameters of the wing are shown in figure 3 a^id the upper 
part of table Vli..L (which follows the form of table 71(a)). The 
chord, (e^^c^), and stiffness matrices have been obtained from the 
given parameters and are shown in the lower part of table v ttt . 

The calculation of the aeroelastlc matrix for the subsonic case 
has been carried out by means of the form of table 7T(b). All but 
three of the steps of the can 5 >utation are shown in table H. nunibered 
in the same order as in table 7T(b). Steps 1, 2, 6, 7, 11, and 12 
ccaistitute matrix multiplications, which are carried out in the order 
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indicated; steps 5 13 constitute matrix additions or outtractlons; 

steps 3 ^ constitute izmltiplications of matrices by constants. 

The aaroelastic matrix is iterated in table X(a) (which follows 
the form of table Vll(a)) to yield a value of ap = -^.208. From this 

value and a value of computed in the same manner for super sonic 

speeds, the divergence dynamic pressure has been calculated by means of 
equation (l6) on the basis of estimated values of the effective lift- 
curve slope. The validation with Mach number of the divergence dynamic 
presGi-ire, the actual dynamic pressure at sea level, and the estimiited 
effective lift-curve slope is shown in figure 4. 


For a value of 


= - 0 . 25 ^ such as would be obtained approximately 


at a Mach number of 1.0, the aerodynamic loading has been calculated for 
the additional angle—of— attack case and the damping— in— roll case in 
table X(b), which follows the form of table VH(b). The values of 


u, 




for the latter case have been calculated from equation (39)* The 
aerodynamic loadings, in addition to those calculated for other 
a/ap values, have been plotted in figure 5 as ratios of the flexible- 
wing loadings to the rigid— wing loadings. The curves have been Integrated 
to yield wing lift, bending-moment , and rolling-moment coefficients as 
well as the aerodynamic center of the wing load, which are shown in 

table X(b) for the case of — = -0.25 and which are plotted against 

aj) 


q/qjj In figure 6. 


The wing lift coefficient Is defined In such a manner that If 
the fuselage lift is known and made dimensionless by dividing by q 
and S the resulting fuselage lift coefficient may be added directly 
to the wing lift coefficient. This definition and the fact that 
figure 5(a) is plotted over the fraction of the wing-alone span 

explains the fact that the area under the curve of figure 5(a) is not 1. 
The aerodynamic center as plotted in figure 6 constitutes the center 
of pressure of only the wing load. To obtain the airplane aerodynamic 
center, the magnitude and center of presaiire of the fuselage load would 
have to be known and taJcen into account. 


Discussion 


Both the aerodynamic and the structural assumqjtions inade in this 
analysis are more realistic than those made in reference 3* The 'device 
employed in this analysis of calculating the air forces for wing sections 
parallel to the direction of fli^t and then transferring them to 
• sections perpendicular to the elastic axis obviates the necessity of 
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replacing 11116 actual wing with one the root and tip of which are 
perpendicular ‘to the elastic axis for the purpose of analysis . Further- 
more , the Inclusion of Falkner's modification (see appendix) In the 
integrating matrices has the effect of rounMng off the load distrihu— 
tion approximately in the manner ohserred at suhsonlc speeds. At 
supersonic speeds the load distrihutions do not go to zero in the manner 
assumed in Fallmer’s modification, hut even at supersonic speeds there 
Is some reduction o£ load at the tip, the totetl magnitude o£ which is 
not far from the reduction obtained hy Fal}mer*s modification. 


Only one aerodynamic assumption is still made, that induction 
effects may he approximated hy an over-all reduction of the strip theory 
loading (rounded off as previously descrlhed) at suhcrltlcal speeds and 
may he neglected at supersonic speeds. The effects of aerodynamic 
Inducticm could he taken into accormt more accurately hy using aero- 
dynamic influence— coefficient matrices instead of the effective lift- 


curve— slope concept and the 


fe] 


and 








matrices used In 


thLis analysis . Available methods of calculating such Influence coeffi- 
cients for vings of arhitrary plan form at suhsonic and supersonic speeds, 
partici^Larly those suitable for wings with large amounts of sweep, are 
either too Inaccurate or too time— cansumlng for practical purposes, 
however. 


"While the analysis of this paper has been performed for wings 
consisting of uncanobered sections, it is directly applicable as well to 
the determlnatian of the additional loading of wings with cambered, 
sections. The loading of such wings due to the section pitching niomeat 
at zero lift may be determined by modifying the analysis somewhat. 

The asBicmptlon of an effective root peirpendlcular to the elastic 
axis made In reference 3 for the purposes of calculating the structural 
response is carried over in this analysis. It is modified, however, 
to the extent that the root is no longer considered to be rigid as in 
reference 3 ^ flexible, both in torsion and bending. It has been 
demonstrated in reference 5 that the deflections of a swept beam may be 
estimated on that assun^jtlon, provided the roob-rotatlon parameters are 
taLown. By assuming the effective root at the Intersection of the elastic 
axis with the side of the fuselage, the root bending due to bending 
moment and root twist due to torque are minimized- The bending due to 
twist axLd twist due to bending are the same regardless of the location of 
the effective root. 

The method of introducing the root rotations into the analysis by 
means of the [l,^ matrix assures that the structural twist in planes 
parallel to the direction of flight is zero at the fuselage. ’Wrani 
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figure 2 it is seen that the local values of the root rotation either 
tend to zero at the root or tend to cancel each other. If the root- 
rotation constants are known, the structural deformations can therefore 
he predicted q^ulte accurately hy the assxnnptlone made. 

The manner In which the equations of equilibrium are solved by 
means of the Integrating matrices accounts for the true chord and 
stiffness variations. It does not necessitate replacement of the 
actual vlng by constant— chord segments with all the flexibility 
concentrated at the ends of the segments, an approach which has been 
used extensively in the work on aeroelastlc problems of strai^t wings. 

A further refinement which obviates the necessity for making any 
structural assumptions other than that of small deflections is the use 
of measured Influence coefficients In the aeroelastlc analysis. Wherever 
such coefficients are available it is, of cotirse, to advantage to use them. 

No explicit account has. been taken In the analysis of the effects 
of the Inertia loading on the structixral deformations and hence the 
aerodynamic loading. On swept wings, In particular, their effects 
may be considerable. For the purposes of this analysis, however, the 
structiiral defoliations due to Inertia loading may be considered part 
of the geometric angle of attack and the rigid-wing geometric angle of 
attack may be modified accordingly. The defoliations due to the inertia 
loading may, IncldenteLlly, be calculated conveniently by means of the 
K 2 , and E^ matrices. 

Some of the general observations made in reference 3 concerning the 
divergence phenomenon are corroborated by the example. As expected of 
a wing with a considerable amount of sweepback, the divergence dynamic 
pressure is negative. Consequently the wing cannot diverge. The 
divergence dynamic pressure is useful as a reference value, however; 
the values of the load distribution and the stability parameters divided 
either by the corresponding rigid-wing values or by the section lift- 
curve elope depend only on the ratio of the actual to the divergence 
dynamic pressure. 

The type of plot shown in figure 4 is therefore quite useful in 
the analysis of aeroelastlc phenomena. As pointed out in reference 3, 
this chart may also be used to estimate the actual divergence dynamic 
pressure where there is a possibility that the wing may diverge. It 
appears that the critical values will occur at either extremity of the 
transonic region. In the transonic region proper the lift-curve slcpe 
usually appears to he lower than at the extremities, so much so that the 
decrease in lift-curve slope even tends to overbalance any forward 
shift in aerodynamic center. 

As would be expected qualitatively, the effect of wing flexibility 
in the case of the example wing is to unload the wing tips due to the 
fact that tl-iey bend up. The lift carried by the wing is therefore less 
than that carried by a rigid wing, the center of pressure being farther 
inboard and the aerodynamic center farther forward. 
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Ihe difference "between the su-pereonlc and subsonic values of the 
loading, the lift and rolliiig-moment coefficients, and the aerodynamic 
center for a given value of • a/a^ is due to the difference In the 

e^ distributions. If the distributions were the same, the subsonic 
■Old supersonic variations of these quantities with a/ajj would coincide 
lesplte the difference in the ej^ values. 

Another Item of possible interest Is the fact that the variations 
of the parameters djj and a^ for the example problem are approximately 

linear (see fig. 7), as would be e3i>ected from, the results of the 
analysis of reference 3* ^Jhe deviations from linearity are most pronounced 
near the points for d = 0 (A = O). Ihey are due to the effects of the 
root rotations, in paarticxilar, the bending due to torsion and torsion due 
to bending; these effects were neglected in the approximate analysis of 
reference 3- Bie points of figure T correspond to the exanqale wing and 
the wings which would be obtained by rotating the example wing to the 
unswept and 37*5° sweptforward positions in such a manner as to keep the 
Cj, cos A (EI) 

parameters — , as well as the chord, stiffness, a.nfl 

®A 

moment-€irm (e^) dlstrlbutlana constant. Points are shown for both the 
subsonio and supersonic variations as well as for the case whai = 0 
over the entire span (ajj = O). The difference between the subsonic and 

supersonic lines is due entirely to the difference in the e^ distribu- 
tion; if the distributions were the same, as woxild be the case if the 
elastic axis were at a constant fraction of the chord, the variations 
would be the same, regardless of the difference in the e^^ values. 

The preseSb analysis is concerned only with wing or tail loads; 
the total loads ar^ obtained by adding the fuselage loads (which may 
be asstaned to be unaffected by flexibility) to the wing oac tall loads 
obtained from the analysis. The amount of load carried by a flexible 
wing and the manner of its distribution can consequently be estimated by 
■the me'fchod presented herein if -the con-taribution of the fuselage is known 
at low dynamic pressures, that is, for the "rigid-firing” case. 

The fuselage has a considerable effect on some of the stability 
parameters as well, al'bhough in -the case of others, such as C 2 , 'the 

effect is negligible. Other effects -that may have to be accounted for 
in calculating s'bability derl'ra.tlves are 'the boundary— layer beha'vior 
and tip suction. The boiindary— layer effect may be accounted for by 
using a section lift— curve slope corrected for boundary— layer effects 
to calculate "the angle— of— attack distribution of the flexible wing at 
the fli^t conditions of Interest and ■then obtaining the lift and drag 
distributions corresponding to that angle-of-attack distribution, 
lateral tip suction may be iBportant on low-aspect -ratio and highly 
swept ■wings. Since it does not affect ■the lift distribution, it may be 
taken into account by calculat in g ■the angle-of-attack distribution of 
the flexible "wlng and estimating the tip suction corresponding to ■fhe 
actual angle of attack at ■fhe tip. 
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In calculating statllity derlTatiTss it Is wall to keep in Tnin<i 
tliat the method presented in this paper is based on a modified strip 
theoi^r, unless aerodyoamlc Influence— coefficient matrices are uaed. Eie 
calculated deriTatires may therefore he samewhat in error, particularly 
if in calculating them the moment of a load dlstrihution has to he 
determined. If there is reason to sxzspeot that the modified strip 
theory is Inadequate for calculating a given derivative, the derivative 
may he calculated for the rigid-wing case hy a more refined method; the 
results calculated hy the method of this paper may then he used to 
correct the accurate rigid-wing value for the effect of structural 
fleiihility. 


COWCinDIlfG EEMAEES 


A method has heen presented for calculating the aerodynamic loading, 
the divergence speed, and certain stahility derivatives of wings and 
tail surfaces of arbitrary plan tona. and stiffness. Provisions have 
been made for using either stiffness curves and root-rotation constants 
or influence coefficients in the structural part of the analysis. Strip 
theory with over-all reduction and rotanding off at the tip to take 
accoimt of aerodynamic Induction have heen used for the aerodynamic part 
of the analysis. Conputing forms, tables of numerical constants required 
in the analysis, and an illustrative example are Included to facilitate 
calculations hy means of the method. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APEERDK 

SUMMAET OF MAIEOZ ALCSBEA PEETINIHT TO TEE AKAIISIS 


For the oonTenience of the reader unfamiliar with matrix terminology, 
a summary of matrix definitions and methods is presented in the following 
sections. For a more oonqplete discussion of matrix methods the reader 
is referred to any. text cn matrices, for instance reference 8. 


Definitions 

A matrix is a rectangular array of numbers, called elements, written 
down in rows and columns . A column matrix consists of a single column, 
a row matrix of a single row. A square matrix has as many rows as it 
has columns. The diagonal of a square matrix from the upper left to the 
lower ri^t is called the principal diagonal . A matrix all the elements 
of which are zero except for those on the principal diagonal is called a 
diagonal matrix. If aU of these elements are unity, the matrix is 
termed the unit matrix. 


Matrix Algebra 

Addition .— Two matrices can be added or subtracted if both hare the 
same number of rows and columns. . The addition or subtraction is carried 
out by adding to or subtracting from each element of the first matrix 
the corresponding element of the second matrix. 

Multiplication by a constant .— A matrix is multiplied by a constant 
by multiplying each element by that constant . 

Matrix multiplication .— Two matrices can be multiplied by each other 
if the second has as many rows as the first ha-s columns. The elements 
of the resulting matrix are obtained by multiplying the elements in the 
corresponding row of the first matrix by those of the corresponding 
column of the second matrix in the following order: The first element 

of the row is multiplied by the first element of the column, the second 
by the second, and' so forth. The sum of the products obtained in this 
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maimer is the Teilue of the element of the produot matrix. Schematically 
this process may he illustrated as follows; 


123456 12345 6 123456 
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Q = aA + hE + cC + dE + eE + f J?' 

It must he eiqjhasized that in multiplying matrices hy each other their 
order is of ln 5 )oid:ance. As the two matrices under consideration are 
written the matrix at the left (the m matrix) is said to he post- 
multiplied hy the other, (the M matrix); or the M matrix may he said 
to he prenniltiplied hy the m matrix, in order to distinguish the manner 
in which they are multiplied. If the two matrices were written in the 
reverse order and then multiplied according to the foregoing instructions, 
that is, if the [m] matrix were postmultiplled hy the [mJ matrix, the 
element of the second row and fourth column oif the product matrix [m] [m] 
would clearly not have the value Q in general, nor would in general 
any other element have the value it would have if the two matrices were 
multiplied in the order shewn. Consequently it is important to observe 
the order in which the matrices are written down in the conputing 
instructions. 

Matrix iteration .— The purpose of iterating a square matrix is to 
determine the column matrix or matrices which, if postmultiplled hy the 
given square matrix, yield the same column matrix except for a canstant 
multiplier. It is the value or values of these multipliers which 
constitute the desired characteristic values of the matrix. 

The iteration is carried out hy assuming a "trial” column, (the 
column shown in table VH is convenient for the purpose of this analysis) 
and premultiplying it hy the given square matrix to yield a "result" 
column. The elements of the result column including the last are 
divided by the last element of the result column, and entered as a 
second trial column. Tlie second trial column is then premultiplied 
hy the square matrix to yield a second result column. The procedure is 
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repeated until the same value (vltliln the desired accuracy) Is obtained 
twice in succession for the last element of the result matrix. Ihe 
reciprocal of this value is the desired (lowest) characteristic value 
of the matrix, that is, the lowest critical value of ajj, in the 
analysis of this paper. 

Another way of estimating a first trial column is to add the 
elements in each row of the A matrix, enter the six sums in the first 
result column and treat them, as if they had been obtained by multiplying 
the A matrix by a first trial column. 


Derivation of the Integrating Matrices 

Althou^ familiarity with the derivation of the integrating 
matrices is not essential to the application of the method of this 
paper, an outline of the derivation is presented because of its general 
interest . 

The integrating inatrices used in this paper are baaed on the same 
concept as Siu5>eon*s rule — replacement of the actual fimctlon which is 
to be Integrated by parabolic segments. If the function y has the 
values respectively, at the equally spaced 

points bhe following relations are seen to be 

true for a second-degree parabola passed throu^ the three known points: 


y = + IC^n+l - yn-l)C3c - + |(7n+l “ 27n + 7n-l)C=^ “ C^l) 

/ 5^^= (-i^yn-l-^(|^)yn + (=J^yn+l ^3) 

^)ra-l * (I * (- ^ (A'‘) 
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f (i - x^)j di = (J- i + (0)7n + (j ^yn+1 ^^5) 

^—1 

j (x - xjy <3^ = ^ yn_l + ^)yn-l + (J ^^n+l 

% 

where 

^ ^ - ^-1 = =4i+l - % 


The different integrations orer the paraholic segments may thus 
he performed hy multiplying the glTen y values by the multiplying 
factors indicated in equations (A2) to (a6). 

Since load distributions at subsonic speeds go to zero with infinite 
slope at the tip and the ordinary second-degree parabola furnishes a 
poor approximation to such a distribution, T. M. Falkner has suggested 
that a curve of the type 


y = Aq + A^(l - + AgCl - x)3/2 


(A8) 


be passed throu^ the last three points of the load— distribution curve 
at the tip (x = l). On the basis of the approximation, relations 
equivalent to equations (A1) to (a 6) may be derived. The multiplying 
factors for the last two segments are then based on these equivalent 
expressions rather than those of equations (A2) to (a6). 

The Integrating factors of equations (A2) to (a 6) may be assembled 

directly into integrating matrices. The Ko matrix, for instance, is 

rx ^ 

set up to perform the integration / y dx. If at the upper limit i = 0.1 

and the ten— point matrix (table IXl(b)) la used, the factors 0.04l67, 
0.06667, -0.00833 may be obtained from equation (aU) since = Oj 

= 0.1, = 0.2, and Ax = 0.1; similarly, if for the same case 

■the integration is extended to x = 0.2 as the upper limit, the 
integrating factors 0.03333, 0.13333, and 0.03333 will be obtained frcmn 
equation (A2). These factors constitute the second and third rows of 
the matrix K^; since the integrations are Independent of the y values 

other than the first bhreo, the other y values are multiplied by zero 
in these two rows, in order to extend the integration to x = 0.3 an 
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integratiaa Is perfoimed up to x = 0.2 and another integration, 

using another paraholic segment, is performed from x= 0.2 to x= 0 . 3 » 
Eor the latter Integration Xq_jj_ = 0.2, x^^ = 0.3, = O.^l-, and 

Ax * 0 . 1 , so that equation (A^v) again yields the factors 0.04 i6y, 

0 . 06667 , and - 0 . 00833 . lEie y Talue at x = 0,2 is therefore assigned 
a multiplying factor of O.O 3333 hy the first integration and a factor 
of 0 . 0^167 hy the second, or a total factor of 0 . 07500 * Ihe resulting 
factors are entered in the fourth row of the matrix. All other 
rows are obtained in a similar manner. 


Ihe matrix is set up to perform the integration 


/ 


y dx 


!Ehe 


values of the last row of the tenr-point matrix (table 1(b)) are 

obtained from Falkner’s equivalent of equation (A3) for the curves 
assumed in reference 5^ with = 0.8, »= 0.9, 3 Cq4.i = 1.0, and 

Ax * 0.1, Only the multiply ing factors for the y values at x « 0.8 
and X = 0.9 are listed, since the y value at x = 1.0 (-the wing tip) 
is assumed to be zero in this analysis, so that its multiplying factor 
is immaterial. Hie values of the last row but one are obtained similarly 
from I’alkner*s equivalent of equation (A2). The values of the row 

for ^ =s 0.7 are obtained by using equation' (A3) in the Interval 
X = 0.6 to X = 0.8 and Falkner*s equivalent of equation (A2) in the 


interval 

obtained 
to 0.8 
interval 


X s 0.8 to 1.0. 


Similarly the row for -3- = 0.6 

®A 


is 


by combining the results of equation (A2) for the interval x 
with FaHaier*s equivalent of equation (A2) for the 
X s 0.8 to 1.0. All other rows are obtained in a similar 


0.6 


manner. 


The matrix is set up to perform the integration 


/>- 


where , x is the variable of integration and x^ the value of x 
the lower limit. In applying the integrating f&tors of equations 
throu^ (a 6) to this integration it must be realized that 


Ax, 

at 

(A2) 


y^(x “ Xo)y dx = (xa - Xq) ^ dx + Jli - XjJj dx 


(A9) 


so that the integrating factors far this Integration would be obtained 
by adding (sj^ ~ Xq) times the factors of equation (A2) or (A 3 ) to 
the factors of equation (A5) or (a 6 ), respectively, depending on the 
limits of the integration. The factors for the different 
Begnents (i = 0.8 to 1 . 0 , 0,6 to 0 . 8 , etc.) are then combined for any 
given row (with its given value of x^) in the manner indicated for 
the El matrix to yield the Z^ matrix. 
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The inatrii sums up the torques outhoard of a given point, 

■while the matrix gives the sum of -fehe moments of forces applied 

outhoard of a given point. Nel'bher requires any Integrations In -the 
sense of equations (A2) to (a 6). For "bhe slx-polnt me'fchod •these ■two 
matr^ices are: 
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It will he noted that ■the moment arms ■which comprise ■the K« matrix are 
fractions of s^, so that the matrix must he multiplied by the length 

In order to yield ac^tual moments as stated In equation (21 ). 
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Table I- Values of Ihe Integrating Matrix CK,] 
(a) S IX -Point Solution 
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(b) Ten -Point Solution 







































































































































Table U- Values of the Integrating Matrix LKtl 
(q) 5ix Point Solution 
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"lable in’- Values of Ihe (ntegnating Matrix CKal 


(a) Six Point Solution 
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( b) Ten Fbint Solution 
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Table IV- Values of Ihe Load Co^^/er5ion Matrix CK«3 
(a) Stx Poinf Solution 
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C b ) Ten Fbint Solution 
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Table Y- Values of Ihe Load Conversion Matrix CK?] 
(a) Six Point Solution 
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Table VI- Form for Computation 0=“ Aeroelastic Matrix 
( a) Win® Parameters 
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table. VI. - CONTINUED 
(b) COMPUTINO INSTRUCTIONS 
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TABLE VL - CONCLUDED 
COMPUTING* INSTRUCTIONS (CONCLUDED) 


SUPERSONIC CASE 


















































Table 31.- Form for solution of Aeroelastic Equation 
fa) Divergence (b) Aerodynamic Loading 
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Table Vm- PARAMETERS of Example Wing 
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TABLEJX.- COMPUTATION OF AEROELASTIC MATRIX OF EXAMPLL WING 
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TABLEX.-SOLimON OF AE.R0ELAST1C EQUATION FOR EMMPLE WING (SUBSONIC CA5£) 
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(b) Aerodynamic Looding 
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(a) Assumed rigid -wing loading. 
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(b) Flexible -wing loading. 


Figure 5.- Load distribution o-f example wing- 









